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There are two methods of preparation of ferrimagnetic spinel. If, during the preparation, an 
external magnetic field as high as 300 Oe is applied upon cooling the material is named field-cooled 
(FC). If the applied field is about lOe the material is zero-field cooled (ZFC). To explore the magnetic 
and thermodynamic properties of these materials we consider two-sublattice spin system, defined 
on the bcc lattice, with spin-s^ operators Sf' at the sublattice A site and spin-s® operators SP at 
the sublattice B site, where s^ > s®. The subtle point is the exchange between sublattice A and 
B spins, which is antiferromanetic. Applying magnetic field along the sublattice A magnetization, 
during preparation of the material, one compensates the Zeeman splitting, due to the exchange, of 
sublattice B electrons. This effectively leads to a decrease of the s® spin. We consider a model with 
s® varying parameter which accounts for the applied, during the preparation, magnetic field. 

It is shown that the model agrees well with the observed magnetization-temperature curves of 
zero field cooled (ZFC) and non-zero field cooled (FC) spinel ferrimagnetic spinel and explains the 
anomalous temperature dependence of the specific heat. 

PACS numbers: 75.50.Gg,71.70.Ej,75.10.Dg,75.10.Lp 


I. INTRODUCTION 

The magnetization-temperature and magnetic suscep¬ 
tibility curves for zero field cooled (ZFC) and non-zero 
field cooled (FC) ferrimagnetic spinel display a notable 
difference below Neel temperature [ll-[I^. The (ZFC) 
curve exhibits a maximum and then a monotonic de¬ 
crease upon cooling from T^v, while the (FC) curve in¬ 
creases steeply, shows a dip near the temperature at 
which the (ZFC) curve has a maximum and finally in¬ 
creases monotonically. There is also difference between 
so-called field-cooled-cooling and field-cooled-warming 
procedures [l3 |. 

The specific heat curves for spinel ferrimagnetics show 
sharp peak at Neel temperature, which indicates fer¬ 
rimagnetic to paramagnetic transition, and tiny peak 
at temperature below Neel’s where the magnetization- 
temperature curve has a maximum iii- 

Although the FC and ZFC spinel ferrimagnetics have 
intensively been studied, their magnetic and thermody¬ 
namic properties have not been understood. There is 
not an effective model which in unified way to explain 
the experimental results. 

In the present paper we consider two-sublattice spin 
system, defined on the bcc lattice, with spin-s^ opera¬ 
tors SA at the sublattice A site and spin-s-® operators 
SP at the sublattice B site, where . The subtle 

point is the exchange between sublattice A and B spins, 
which is antiferromanetic. Applying magnetic field along 
the sublattice A magnetization, during preparation of the 
material, one compensates the Zeeman splitting of sub¬ 
lattice B electrons due to the exchange. This effectively 
leads to decrease of the s® spin. One can obtain an intu¬ 
ition for this from spin-fermion model of spinel ferrimag¬ 
netic with spin-1/2 itinerant electrons at the sublattice B 
site and spin-s localized electrons at the sublattice A site. 
An applied, along the magnetization of the localized elec¬ 


trons, external magnetic field compensates the Zeeman 
splitting due to the spin-fermion exchange. Integrating 
out the fermions we obtain an effective spin model with 
effective spin which depends on the external magnetic 
field (see Appendix A). 

We consider a model with s® varying parameter which 
accounts for the applied, during the preparation, mag¬ 
netic field. First we present the method of calcula¬ 
tion exploring ZFC spinel with = 1.5 and = 1. 
We obtained that the system has two phases. At low 
temperature (0,r*) the magnetic orders of the A and 
B spins contribute to the magnetization of the sys¬ 
tem, while at the high temperature (T*,Tn), the mag¬ 
netic order of the sublattice B with smaller spin and 
with a weaker intra-sublattice exchange is suppressed by 
magnon fluctuations. Only the sublattice A spins, with 
stronger intra-sublattice exchange, have non-zero spon¬ 
taneous magnetization. There is no additional symmetry 
breaking, and the Goldstone boson has a ferromagnetic 
dispersion in both phases, but partial-order transition 
demonstrates itself through tiny peak of specific heat as 
a function of temperature. 

Partial order is well known phenomenon and has been 
subject to extensive studies. Frustrated antiferromag¬ 
netic systems has been studied by means of Green func¬ 
tion formalism. Partial order and anomalous temper¬ 
ature dependence of specific heat have been predicted 
[l5| . Experimentally the partial order has been observed 
in Gd 2 Ti 207 [l^. Monte Carlo method has been uti¬ 
lized to study the nature of partial order in Ising model 
on kagome lattice [a. There are exact results for the 
partially ordered systems which precede the above stud¬ 
ies [iTHl^ . The advantage of the present method of cal¬ 
culation is that it permits to consider the sublattice B 
spin as a varying parameter s'® < 1. We calculate the 
magnetization-temperature curves and the specific heat 
as a function of temperature for s® = 0.7, s® = 0.4, 
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= 0.2 and fixed = 1.5, thus accounting for the in¬ 
creasing of the applied, during preparation of FC spinel, 
magnetic field. 


II. METHOD OF CALCULATION 

The Hamiltonian of the ZFC ferrimagnetic spinel is 

H = -j^ E s/ - J"" Y. s? • sf 

+ jJ2sf-sf (1) 

(d> 

where the sums are over all sites of a body-centered cu¬ 
bic (bcc) lattice: {i,j) denotes the sum over the nearest 
neighbors, <C i,j denotes the sum over the sites of 
the A sublattice, <C i,j denotes the sum over the 
sites of the B sublattice. The first two terms describe 
the ferromagnetic Heisenberg intra-sublattice exchange 
> 0, > 0, while the third term describes the inter¬ 

sublattice exchange which is antiferromagnetic J > 0. 

To study a theory with the Hamiltonian Eq.(IT|) it 
is convenient to introduce Holstein-Primakoff represen¬ 
tation for the spin operators Sf(a+,a) and Sf(6+,6). 
Rewriting the effective Hamiltonian in terms of the Bose 
operators (a+,a, 5+,6) we keep only the quadratic and 
quartic terms. The next step is to represent the Hamil¬ 
tonian in the Hartree-Fock approximation: 

H « Hhf = Ha + Hq (2) 

with 

Ha = -if+ QNJ^{s^f{u^-if 

+ mjs^s^{u-lf, (3) 

and 

Hq=Y - Ik {atK+ bkak)] , 

k^Br 

( 4 ) 

where N = is the number of sites on a sublat¬ 

tice. The two equivalent sublattices A and B of the bcc 
lattice are simple cubic lattices. The wave vector k runs 
over the reduced first Brillouin zone Br of a bcc lattice 
which is the first Brillouin zone of a simple cubic lattice. 
The dispersions are given by equalities 

el = (3 — cos kx — cos ky — cos kf + 8s^ Ju 

e\. = 4s^ (3 — cos kx — cos ky — cos kf + Ju 

o T r~A B kx ky kz - , 

7 fe = 8J uvs^ cos—coscos— (5) 

The equations show that Hartree-Fock parameters 
{u^, , u) renormalize the intra and inter-sublattice ex¬ 

change constants {J^, J^, J) respectively. 


To diagonalize the Hamiltonian one introduces new 
Bose fields ak, a'^, Pk, by means of the transforma¬ 
tion 

ak = Ukak + Vk Pk ='^kOi'l + Vk Pk 

( 6 ) 

bk = UkPk + Vka'^ h'l = UkPk + ^kotk, 

where the coefficients of the transformation Uk and Vk are 
real functions of the wave vector k (lB7l) . The transformed 
Hamiltonian adopts the form 

H,= Y + ^k + Ei) , (7) 

fceSr 

with new dispersions 



and vacuum energy 



For positive values of the Hartree-Fock parameters and 
all values of fc € Br, the dispersions are nonnegative 

> 0, > 0. When the ak boson is 

the long-range (magnon) excitation in the system with 
E'^ oc p/c^, near the zero wavevector, while the Pk bo¬ 
son is a gapped excitation, with gap proportional to the 
inter-sublattice exchange constant Eq = 8Ju{s^ — s^) . 

The free energy of a system with Hamiltonian Hhf 
equations and @ is 

E = 6NJ^{s^f{u^-If+6NJ^{s^f{u^-If 

+ 8NJs^s^{u-lf + j^J2^°^ (10) 

^ kGBr 

+ mY (l - + In (l - , 

^ fceSr 

where P = l/T is the inverse temperature. Then, the 
system of equations for the Hartree-Fock parameters is 

dF/du^ = ll, dE/du^ = ll, dFldu = ll. ( 11 ) 

The Hartree-Fock parameters are positive functions 
of T/J , solution of the system of equations (fTTI) (see 
Eas. (IB8|) l. Utilizing these functions, one can calculate 
the spontaneous magnetization on the two sublattices 

= < Sf > j is from sublattice A 

( 12 ) 

= < S 2 j > j is from sublattice B 
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and M = , the spontaneous magnetization of 

the system. In terms of the Bose functions of the a and 
/3 excitations they adopt the form 




^ k^B-r 

^ [^Irik + ul nl + vl 


(13) 


where Uk and Vk are functions of the wavevector fc, coef¬ 
ficients in the transformation ([5]). 

The magnon excitation - is a complicated mixture of 
the transversal fluctuations of the A and B spins ([5]). As 
a result the magnons’ fluctuations suppress in a different 
way the magnetization on sublattices A and B. Quantita¬ 
tively this depends on the coefficients Uk and Vk- At char¬ 
acteristic temperature T* spontaneous magnetization on 
sublattice B becomes equal to zero, while spontaneous 
magnetization on sublattice A is still nonzero (see FigllJ. 

To study the magnetic properties of the system above 
T* we first consider the paramagnetic phase. To this end 
we make use of the Takahashi modified spin-wave theory 
and introduce two parameters and A® to enforce 
the magnetization on the two sublattices to be equal to 
zero. The new Hamiltonian is obtained from the old one 
equation o by adding two new terms: 


H = H - + Y (14) 

ieA ies 


In momentum space the new Hamiltonian H adopts the 
form Eq.(|3]) with new dispersions 

=el + A^, el = el + X^. (15) 

Utilizing the same transformation (|6]) one obtains the 
Hamiltonian H in diagonal form 0 with dispersions E°‘ 
and obtained from Eqs.lH]) replacing e% and el with 
el and el respectively. It is convenient to represent the 
parameters A^ and A® in the form 

A^ = 6Jus^{^"^ — 1), A'® = — 1). (16) 

The dispersions e^ and el are positive for all values of the 
wavevector k, if the parameters and are positive. 

The dispersions E^ and E^ are well defined if 

/iV > 1- (17) 

The /3fc excitation is gapped {E^ > 0) for all values of 
parameters /i^ and which satisfy equation (1171) . The 
a excitation is gapped if > 1, but in the particular 

case 

= 1 (18) 

Eq = 0, and near the zero wavevector E^ « There¬ 
for, in the particular case Eq. (flSl) ak boson is the long- 
range excitation (magnon) in the system. 



FIG. 1: (Color online) Spontaneous magnetization M^, 
and as a function of the temperature in units of 

exchange constant T/J for ferrimagnetic spinel on bcc lattice 
with = 1.5, s® = 1, B/J = 2 and /J = 0.002. The 
partial order transition temperature T* is the temperature 
above which = 0. 


Above Neel temperature we introduced the parameters 
/i^ and fj,^{X"^ , X^) to enforce the sublattice A and B 
spontaneous magnetizations to be equal to zero. We find 
out these and Hartree-Fock parameters, as functions of 
temperature, solving the system of five equations, equa¬ 
tions m and the equations = 0, where the 

spontaneous magnetization has the same representation 
as equations (US but with coefficients Uk, Vk, and disper¬ 
sions E^ in the expressions for the Bose functions. 
The numerical calculations show that above Neel tem¬ 
perature > 1. When the temperature decreases 

the product decreases, remaining larger than one. 

The temperature at which the product becomes equal to 
one = 1) is the Neel temperature Fig. ([6]). 

Below T/v, the spectrum contains long-range (magnon) 
excitations. Thereupon, = 1 and one can use the 

representation = l//i^. Then, and Hartree-Fock 
parameters are solution of a system of four equations, 
equations m and the equation = 0. 

We utilize the obtained functions ^"^(T),/i'®(T) u^{T), 
u^{T), u{T) (see Appendix C) to calculate the sponta¬ 
neous magnetization as a function of the temperature. 
For a system with = 1.5, = 1, J^/J = 2 and 

J^/J = 0.002 the functions M^lT/J), M^lT/J) and 
M^{T/J) + M^{T/J) are depicted in figure ([T]). The 
partial order transition temperature T* is the tempera¬ 
ture above which = 0. 

The customary formula for the entropy of a Bose sys¬ 
tem with Hamiltonian ([7]) is 

5 = ^ ^ [(1 -b uf ) ln(l -b 4) - 4 In 4] , (19) 

k,6 

where 6 stays for a and (3. The dispersions E^ and E^ 
([51) are used to define the Bose functions and nf below 
T *, dispersions E^ and E^ with = 1 are used for 
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FIG. 2: (Color online) Specific heat C vs temperature, in 
units of the exchange constant T/J, for j J = 2, / J = 
0.002, = 1.5 and s® = 1. The high temperature (green) 

vertical line marks Neel Tjv temperature, while the low tem¬ 
perature (red) line marks partial order transition temperature 

partial order phase T* < T < Tat, and with > 1 

for paramagnetic phase above Neel temperature. With 
entropy, as a function of temperature in mind, one can 
calculate the contribution of magnons to the specific heat: 

C = T§ (20) 

The resultant curve C{T/J), for a system with the same 
parameters, as above, is depicted in figure dH). 

The figures o and ([2]) show that the present method of 
calculation describes correctly the features of the system, 
partial order transition and anomalous behavior of the 
specific heat at the temperature of this transition T*. 

III. FC SPINEL FERRIMAGNETICS 

FC spinel ferrimagnetics are prepared applying mag¬ 
netic field along the sublattice A magnetization upon 
cooling the material. This compensates the Zeeman split¬ 
ting of sublattice B electrons which effectively leads to 
decrease of the spin. The Hamiltonian of the FC 
spinel ferrimagnet is given by Eq. © with a varying 
parameter. The advantage of the method of calculation, 
presented above, is that we can use the same systems of 
equations for the three phases, 0<r < T*, T* < T < 
Ta/, T > Tjv, and different values of the parameter . 
We calculate the magnetization-temperature curves and 
the specific heat as a function of temperature for = 1.5 
and three different values of {s^ = 0.7, = 0.4 and 

= 0.2), thus accounting for the increasing of the ap¬ 
plied, during preparation of FC spinel, magnetic field. 

The resultant magnetization-temperature curves are 
depicted in figure ®. The curve ”a” corresponds to 
ZFC spinel (see hgure ©)• The increasing of the ap¬ 
plied, during the preparation, magnetic field is mod¬ 
eled by decreasing of . The curves ”b”, ”c” and ”d” 



FIG. 3: Color online) Spontaneous magnetization 
as a function of the temperature in units of the exchange 
constant T/J for ferrimagnetic spinel on bcc lattice with 
J'^/J = 2, J® = 0.002, = 1.5 and :a)(black) s® = 1, 

b)(red) s® = 0.7, c)(blue) s® = 0.4 and d)(green) s® = 0.2. 
The dash (magenta) line shows the spontaneous magnetiza¬ 
tion of sublattice A for = 1.5 and s® = 0.2. 



FIG. 4: (Golor online) Specific heat C vs temperature, in 
units of the exchange constant T/J, for = 1.5 and different 
values of the effective spin s^: black curve s® = 1, red curve 
s® = 0.7 and blue curve s® = 0.2 . 


are magnetization-temperature curves for FC spinel fer¬ 
rimagnetics with S'® = 0.7, s® = 0.4 and s® = 0.2. The 
curve ”c” agrees well with the observed magnetization- 
temperature curve of FC spinel MnV20i with an exter¬ 
nal magnetic field as high as 300 Oe applied upon cool¬ 
ing the material ii- The dash line shows the spon¬ 
taneous magnetization of sublattice A for s"^ = 1.5 
and s® = 0.2. Comparing with line ”d”, which shows 
the spontaneous magnetization + M® of the same 
system one can conclude that contribution of the sublat¬ 
tice B magnetization is small and Zeeman splitting, of 
sublattice B electrons, is approximately compensated. 

The dependence of specific heat on temperature for 
different values of s® is shown in figure ([1]). The curve 
with higher maximum corresponds to ZFC spinel (see fig¬ 
ure ©) With increasing the applied, during the prepara¬ 
tion, magnetic field (with decreasing s®) the anomalous 
temperature behavior of specific heat, at temperature of 
partial order transition, remains well defined but the pick 
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decreases. The anomalous pick is well observed experi¬ 
mentally d, 0,0 , and one can use it to determine the T* 
temperature of the partial order transition. 

IV. SUMMARY 

In this paper is pointed out that a model which agrees 
well with the observed magnetic and thermodynamic 
properties of ZFC and FC spinel ferrimagnets is a two- 
sublattice Heisenberg model with spin-s^ operators at 
the sublattice A site, spin-s^ operators at the sublat¬ 
tice B site, ferromagnetic intra-sublattice exchange and 
antiferromagnetic inter-sublattice exchange. The ap¬ 
plied magnetic field, along the sublattice A magnetiza¬ 
tion upon cooling the material, is accounted for varying 
only one parameter . 

The subtle point is that MnVzOi spinel has an obvious 
anomalous magnetic behavior, but the appearance of this 
anomaly is attributed to the strong orbital-spin coupling 
[3 which is not discussed in the model (H]). The spinel 
MnV 204 is a two-sublattice ferrimagnet, with site A oc¬ 
cupied by the M'n?^ ion, which is in the high-spin 
configuration with quenched orbital angular momentum, 
which can be regarded as a s = 5/2 spin. The B site is 
occupied by the ion, which takes the high-spin 
configuration in the triply degenerate t 2 g orbital and has 
orbital degrees of freedom. Because of the strong spin- 
orbital interaction it is convenient to consider jj coupling 
with and J® = . The sublattice A 

total angular momentum is = 5/2, while the 

sublattice B total angular momentum is + s^, 

with = 3, and s'® = 1 [I|. Then the g-factor for 
the sublattice A is = 2, and for the sublattice B 
= |. The sublattice A magnetic order is antiparallel 
to the sublattice B one and the saturated magnetization 
is (T = 2| — |4 = 0, in agreement with the experimen¬ 
tal finding for ZFC spinel that the magnetization goes to 
zero when the temperature approaches zero. The Hamil¬ 
tonian of the system is 

H = Y. E Jf-Jf 

<ij»A 

+ ( 21 ) 

(b) 

The first two terms describe the ferromagnetic Heisen¬ 
berg intra-sublattice exchange > 0 ,/c® > 0, while the 
third term describes the inter-sublattice exchange which 
is antiferromagnetic k > 0. The operators 3^ and J® sat¬ 
isfy the SU{2) algebra, therefor we can use the Holstein- 
Primakoff representation of the total angular momen¬ 
tum vectors 3^{a'^,aj) and J® (6^, bj), where a^, aj and 
6^, bj are Bose fields. Farther on, we repeat the calcula¬ 
tions from sections H and HI. The only difference is the 
expression for the magnetization + g® M®. The 

jj model shows that the anomalous magnetic and ther¬ 
modynamic behavior of MnV 204 is a consequence of the 


ferrimagnetic nature of the system. 

Appendix A: Spin-fermion model of spinel 
ferrimagnet 

The Hamiltonian of the spin-fermion model of ferri¬ 
magnetic spinel defined on a body centered cubic lattice 
is 

H = - t Y + h.C.) - /i^ Ui 



i^B 


j®^ E 

sp-sf+ jEs^- 

Sf (Al) 

<^ij'»B 

{ij) 



t-sf-HYsA 

-hYsA, 


i^A 

ieB 


where S'/® = 5 X) ^ with the Pauli matrices 

<7 < 7 ' 

,T^), is the spin of the itinerant electrons at the 
sublattice B site , SA is the spin of the localized electrons 
at the sublattice A site, g, is the chemical potential, and 
rii = cf^Cia- The sums are over all sites of a body cen¬ 
tered cubic lattice, (i, j) denotes the sum over the nearest 
neighbors, while <C ij and ij are sums over 
all sites of sublattice A and B respectively. The Heisen¬ 
berg term {J^ > 0) describes ferromagnetic Heisenberg 
exchange between localized electrons and J > 0 is the an¬ 
tiferromagnetic exchange constant between localized and 
itinerant electrons. H > Q is the Zeeman splitting en¬ 
ergy due to the external magnetic field (magnetic field 
in units of energy). We represent the Fermi operators, 
the spin of the itinerant electrons and the density oper¬ 
ators riio- in terms of the Schwinger bosons 
and slave fermions {hi,hf ,di,df). The Bose fields are 
doublets (tr = 1,2) without charge, while fermions are 
spinless with charges 1 (di) and -1 (hi): 

Cif = hf Cii = hf (pi2 - 

ni = I- Kl hi + d+ di, s/ = ^ E ’ 

<7 < 7 ' 

= df di (A2) 

+ d'ldi + hf hi = l (A3) 

To solve the constraint fEa lA3l) . one makes a change of 
variables, introducing Bose doublets Qa and [2l[ 

Cicr — g^i(7 (1 h^ hi d^ di^ , 

c+ = ^tA^-dth^-dtd,y\ (A4) 

where the new fields satisfy the constraint Cil/Cio- = 1- In 
terms of the new fields the spin vectors of the itinerant 
electrons have the form 

sy = i E [1 - Ah^- dfd,] (A5) 

( 7 ( 7 ' 
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When, in the ground state, the lattice site is empty, the 
operator identity hfhi = 1 is true. When the lattice 
site is doubly occupied, df di = 1. Hence, when the lat¬ 
tice site is empty or doubly occupied the spin on this 
site is zero. When the lattice site is neither empty nor 
doubly occupied {hf hi = df di = 0), the spin equals 
Si = l/2ni, where the unit vector 

(n? = l) (A6) 

gg' 

identifies the local orientation of the spin of the itinerant 
electron. 

The part of Hamiltonian Eq. m with itinerant 
fermions can be rewritten in terms of Bose fields Ea. (IA4l) 
and slave fermions 

= -t Md^-hp,)CtaCJ. 

+ {dpi - d+h+) (Ca02 - O2O1) + h.c.] 

X (1 - h+h, - d+d,) ^ (1 - h+hj - d+d,) ^ 

— E rriimj rii • nj 

(1 - h+h, + d+di) (A7) 

ieB ieB 

les ^ 

where 

mi = ^[l- hfh^- df di\. (A8) 

The mixed, spin-fermion term, adopts the form 

= (A9) 

hi) 

The Hamiltonian (lAll) can be rewritten in the form 

H = H^ + H^ + (AlO) 

where is the contribution of sublattice A spins SA. 

An important advantage of working with Schwinger 
bosons and slave fermions is the fact that Hubbard term 
is in a diagonal form. The fermion-fermion and fermion- 
boson interactions are included in the hopping term. 
One treats them as a perturbation. To proceed we ap¬ 
proximate the hopping term of the Hamiltonian Eq. (IA7I) 

setting {1 — hfhi —d'idi)'^ ^ 1 and keeping only the 
quadratic, with respect to fermions, terms. This means 
that the averaging in the subspace of the fermions is per¬ 
formed in one fermion-loop approximation. Further, we 
represent the resulting and Hamiltonian as a sum 
of two terms 

Hf = Hi+ H(i, 


where 

hI = - t E (didi - hih^ + h.c.) +U'^didi 

<.ii:$>B i&B 

- ^JiY,{l-htK + dtdi) (A12) 

i&B 

ieB 

is the Hamiltonian of the free d and h fermions, and 

hL= - t E [(did.-hih;){CtQ^-l){AU) 

+ {dihi - d+hi) (C.iC ,2 - foOl) + h.c.] 

— ''PI rriimj rii • nj 

<^ij'»B 

is the Hamiltonian of boson-fermion interaction. 

The ground state of the system, without accounting for 
the spin fluctuations, is determined by the free-fermion 
Hamiltonian ho and is labeled by the density of electrons 

n = 1— < hf hi > + < df di > (A14) 

(see equation (|A2I) ') and the ’’effective spin” of the sub¬ 
lattice B electron 

= \ (1- < hth^ > - < d+d* >) . (A15) 

At half-filling 

< hi hi >=< didi> . (A16) 

To solve this equation one sets the chemical potential 
/i = U/2. Utilizing this representation of fi we calculate 
the effective spin as a function of applied magnetic 
field h for parameters At/U = 2.2 and As^J/U = A. The 
result is depicted in figure ([5]). 

Let us introduce the vector, 

Ali = E = pf. (A17) 

gg' 

Then, the spin-vector of itinerant electrons Ea. (IA5ll can 
be written in the form 

(1 - hi h, - di di) , (A 18 ) 

where the vector identifies the local orientation of the 
spin of the sublattice B itinerant electrons. 

The Hamiltonian is quadratic with respect to the 
fermions di , di and hi,hi, and one can average in the 
subspace of these fermions (to integrate them out in the 
path integral approach). As a result, accounting for the 
definition of (see Ea lAlSI) . one obtains and 

an effective model for spin vectors with Hamiltonian 

<.i3:>B 


(All) 


(A19) 
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FIG. 5: (Color online) Sublattice B effective spin s® as a 
function of applied magnetic field h for parameters At/U = 2.2 
and As^J/U = 4. 


where the effective exchange constant is calculated 
in the one loop approximation. The term with ex¬ 
change constant in equation (IA7D and the mixed, spin- 
fermion term (|A9I) . adopt the form 

(A20) 

and 

(A21) 

(ij) 


respectively. Collecting all terms one obtains the Hamil¬ 
tonian of the FC ferrimagnetic spinel 

H = -J^ Y • s/ - Y S? • Sf 
+jY^t-sf, 

<b> 


where and s® depending on the field 

applied during preparation Fig.([5|). After the process of 
preparation the external magnetic field is set equal to 
zero. For systems with = 1, discussed in the pa¬ 
per, one has to consider two-band model for sublattice 
B fermions. The result remains the same even for this 
more complicate system: the effective spin decreases 
when the applied, under the preparation, magnetic field 
increases. 


Appendix B: Hartree-Fock approximation 

Let us consider a theory with Hamiltonian (IT|). We 
introduce Holstein-Primakoff representation for the spin 


operators 

sf+ = Sf + iSf = - a+a, a, 

Sf- = Sf'-iSf=a*^2s^-a*n, (Bl) 

when the sites j are from sublattice A and 

SY = sf^ + ^Sf^ = - 6 + ^2sB - bp, 

Sf- = Sf^ - zSp = -^2sB - bp, b, (B2) 

sr= 

when the sites j are from sublattice B. The operators 
a^, aj and 6^, hj satisfy the Bose commutation relations. 
In terms of the Bose operators and keeping only the 
quadratic and quartic terms, the effective Hamiltonian 
Eq.dl]) adopts the form 

H = H2 + H4 (B3) 


where 

H2 = Y, - a^cLi - afaj) 

+ Y. + b+bj - bp, - b+bj) (B4) 

+ JYI {af 6+ -I- aibj) 

<b> 






+ Y [bfbpbi-bjf+ {b+-bp^hbj] 




<y> 
7b 


{a,b+ bjbj + a+ b+ b+ bj) (B5) 


-by — {af ttittibj + af af ap ) — Ap a^b^ bj 


The next step is to represent the Hamiltonian in the 
Hartree-Fock approximation ismh and 


= ■£ 

{api + apj - api - apj) 


+ Y 

{bp, + bp, - bp, - bp,) {B6) 

<.ij:^B 

+ JuYls'^b+bj 

+ s^api - V ipb'^ + aib,) 

{b> 



It is convenient to rewrite the Hamiltonian in momen¬ 
tum space representation Eq.(|4]) The two equivalent sub¬ 
lattices A and B of the bcc lattice are simple cubic lat¬ 
tices. The wave vector k runs over the reduced first Bril- 
louin zone Br of a bcc lattice which is the first Brillouin 














zone of a simple cubic lattice. The dispersions are given 
by equalities Eqs. m 

To diagonalize the Hamiltonian one introduces new 
Bose fields /3fc, by means of the transforma¬ 

tion (|6]) where the coefficients of the transformation Uk 
and Vk are real function of the wave vector k 


Uk 


Vk 


if 



sj{el + elY - 472 ) 


(B7) 


sign{-fk) 



The transformed Hamiltonian adopts the form Eqs. o 
with dispersions Eqs.® The free energy of a system with 
Hamiltonian Hhf (di) is Eq. (na The system of equa¬ 
tions for the Hartree-Fock parameters Eqs. dm , in terms 
of the Bose functions of the a and /? excitations, adopt 
the form 


u = 


,_ 

3sA n ^ 
k£Br 

1 - 


3s^ N 


2 a I 2 8 , 2 

^k ^k ^k ^k 


2a, 20, 2 

Vk Uk + u^nl + Wfc 


feGSr 


= 1 - 




k£Br 


+ ^ < + Wfe nf + 


- SJu[l + nl+n^0j 


^cos ^ cos ^ cos 


{el + elf - 


21 


where and are the Bose functions of a and /3 ex¬ 
citations. 


Appendix C: Takahashi modified spin-wave theory 


To study the paramagnetic phase of the system we 
make use of the Takahashi modified spin-wave theory 
and introduce two parameters and to enforce the 
magnetization on the two sublattices to be equal to zero. 
It is convenient to represent the parameters A"^ and A^ in 
the form Eas. dTBl) . We find out these and Hartree-Fock 
parameters, as functions of temperature, solving the sys¬ 
tem of five equations, equations (IB8I) and the equations 
= 0. The numerical calculations show that 


above Neel temperature > 1 (blue line in Fig.®). 

When the temperature decreases the product de¬ 

creases, remaining larger than one. The temperature at 
which the product becomes equal to one = 1) is 



FIG. 6: (Color online) Parameters p® and as a 

function of the temperature, in units of the exchange con¬ 
stant r/J, for ferrimagnet with sublattice A spin = 1.5 
and sublattice B spin s® = 1. The high temperature vertical 
(cyan) line corresponds to the Neel temperature of ferrimag¬ 
net to paramagnet transition. The low temperature vertical 
(green) line corresponds to partial order transition tempera¬ 
ture. 


the Neel temperature, marked by vertical cyan line in 
Fig.®. 

At low temperature = 1 and = 1 (A"^ = A^ = 
0). The Hartree-Fock parameters are positive functions 
of T/J , solution of the system of equations (IB8I) . Uti¬ 
lizing these functions, one can calculate the spontaneous 
magnetization on the two sublattices at low temperature. 
At characteristic temperature T* spontaneous magneti¬ 
zation on sublattice B becomes equal to zero, while spon¬ 
taneous magnetization on sublattice A is still nonzero 
(see Fig.l in the paper). 

Above T* the spectrum contains long-range (magnon) 
excitations. Thereupon, = 1 and one can use the 

representation ■ Then, and Hartree-Fock 

parameters are solution of a system of four equations, 
equations (IB8I) and the equation = 0. 

The functions fi^(T/J), g^(T/J) and 

J)g.^ {T/J) are depicted in figure ® for a 
system with = 1.5, S'® = 1, J^/J = 2 and 
JB/J = 0.002. 

For the same system the Hartree-Fock parameters 
u^{T/J), u^{T/ J) u{T/ J) are depicted in figure ®. 
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